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Kinetic Boltzmann equations are used to model the ionization and expansion dynamics ol xenon 
clusters irradiated with short intense VUV pulses. This unified model includes predominant in- 
teractions that contribute to the cluster dynamics induced by this radiation. The dependence of 
the evolution dynamics on cluster size, N at0 ms = 20 — 90000, and pulse fiuence, F — 0.05 — 1.5 
J/cm 2 , corresponding to intensities in the range, 10 12 — 10 14 W/cm 2 and irradiation times, < 50 
fs, is investigated. The predictions obtained with our model are found to be in good agreement 
with the experimental data. We find that during the exposure the cluster forms a shell structure 
consisting of a positively charged outer shell and a core of net charge equal to zero. The width 
of these shells depends on the cluster size. The charged outer shell is large within small clusters 
{Natoms = 20,70), and its Coulomb explosion drives the expansion of these clusters. Within the 
large clusters (N a toms = 2500, 90000) the neutral core is large, and after the Coulomb explosion of 
the outer shell it expands hydrodynamically. Highly charged ions within the core recombine effi- 
ciently with electrons. As a result, we observe a large fraction of neutral atoms created within the 
core, its magnitude depending on the cluster size. 

PACS numbers: 41.60.Cr, 52.50.Jm, 52.30.-q, 52.65.-y 



Atomic clusters are excellent objects to test the dy- 
namics within samples irradiated with radiation from 
short wavelength free-electron-lasers (FELs) P, H, ||. 
Their physical properties put them on the border be- 
tween the solid state and the gas phase. Cluster stud- 
ies are important for planned experiments with FELs in 
solid state physics, materials science and for studies of 
the extreme states of matter Q ■ Accurate predictions on 
the ionization, thermalization and expansion timescales 
within irradiated samples that can be obtained with clus- 
ter experiments are also needed for exploring the lim- 
its of experiments on single particle diffraction imaging 

During the first cluster experiments performed at the 
free-electron-laser facility FLASH at DESY with VUV 
photons of energy, E = 12.7 eV, and power densities 
up to a few 10 13 W/cm 2 [ll[ highly charged Xc ions 
(up to +8) of high kinetic energies were detected. This 
unexpectedly strong energy absorption could not be ex- 
plained using the standard approaches [l]| 0, Q3] ■ More 
specifically, the energy absorbed was almost an order of 
magnitude larger than that predicted with classical ab- 
sorption models, and the ion charge states created were 
much higher than those observed during the irradiation of 
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isolated atoms in similar conditions. This indicated that 
at such radiation wavelengths some processes specific to 
many-body systems were responsible for an enhanced en- 
ergy absorption. 

The physics underlying the dynamics within the ir- 
radiated clusters is complex. Several interesting theo- 
retical models have been proposed in order to describe 
the evolution of clusters exposed to intense VUV pulses 
[H [H, El GJ, HUH. The most explored ones are: 
i) heating of quasi-free electrons due to enhanced inverse 
bremsstrahlung (IB) [IE], [l6j] , ii) enhanced photoioniza- 
tion within the sample due to lowering of interatomic 
potential barriers [17], [HI > and iii) heating due to many- 
body recombination processes [2fJ] . Each of these various 
approaches lead to the significantly enhanced energy ab- 
sorption in agreement with the experimental data [ll[. 
However, one would expect that if all proposed enhance- 
ment factors are included within one model, it would 
probably lead to absorption rates much higher than those 
experimentally observed. 

In order to evaluate the contribution of various pro- 
cesses to the ionization dynamics, we have constructed a 
unified model [13] , [23[ based on kinetic equations includ- 
ing the following predominant interactions: photoion- 
ization, collisional ionization, clastic scattering of elec- 
trons on ions, inverse bremsstrahlung heating, electro- 
static interactions between charges and with laser field, 
shifts of energy levels within atomic potentials due to 
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the plasma environment, and shielded electron-electron 
interactions. Using a non-equilibrium Boltzmann solver, 
we followed the evolution of the Xe clusters of various 
size (N a toms = 20 — 90000) irradiated with rectangular 
VUV pulses of intensity 10 12 — 10 14 W/cm 2 and duration, 
< 50 fs. We found that all physical mechanisms that were 
included into the model contributed to the ionization dy- 
namics but with different weights. The total ionization 
rate within the sample was affected most by the inverse 
bremsstrahlung heating rate applied. Within the the- 
oretical framework defined above we estimated that: i) 
many-body effects (many-body recombination) 20] could 
contribute only for clusters irradiated at low pulse flu- 
ences, ii) the plasma environment effects estimated with 
electron screened atomic potentials were small. This was 
in contrast to the estimates of Ref. [fTf performed with 
unscreened atomic potentials. Our study extended the 
treatment of Ref. [161 ] , where neither the spatial inhomo- 
gcneity of the clusters was treated nor the gross move- 
ment of electrons. It turned out that both effects sig- 
nificantly contribute to the cluster dynamics: the struc- 
ture of charge created within the cluster was found to be 
strongly inhomogeneous (24[. This inhomogeneity was 
induced by the dynamics of electrons. 

The main aim of this letter is to obtain a complete de- 
scription of ionization and expansion dynamics. There- 
fore we have included into the present analysis the evalua- 
tion of the effect of the three-body-recombination process 
(the non-equilibrium recombination rates are included 
into kinetic equations), and the expansion of the clus- 
ter until ions start to leave the simulation box (up to ~ 2 
ps). 

Dynamics within irradiated clusters depends on the 
pulse flucnce and cluster size. We performed simulations 
of irradiated clusters at the parameter values correspond- 
ing to those at the first FLASH cluster ex peri ment at 
DESY at 100 nm radiation wavelength [Hill- These 
predictions will now be compared to the experimental 
data. 

The simulations were performed for clusters exposed 
to single rectangular VUV pulses of a fixed fluence, but 
of various intensities and pulse durations. The pulse 
intensity was 10 12 — 10 14 W/cm 2 and the pulse length 
was, At < 50 fs. The predictions obtained from differ- 
ent events were then averaged over the number of events. 
This scheme followed the experimental data analysis: ex- 
perimental data were obtained after averaging the single 
shot data obtained with FEL pulses of various temporal 
shapes but of a fixed radiation flux. In particular, the 
experimentally estimated histograms of ion charge and 
predictions on the average kinetic energy per ion were 
obtained from the averaged time-of-flight (TOF) spec- 
tra. The TOF detector could record charged particles 
(ions) only. The experimental ion intensities were ob- 
tained by integrating and averaging the TOF signal over 
subsequent FEL pulses. Some of those intensities were 
corrected for the relative geometric acceptances of the 
TOF detector or the MCP detector efficiencies for differ- 



ent charge states. The maximal experimental uncertainty 
of the fluence estimation is given by factor of 2.5. We ob- 
tain predictions for two limiting cases: i) a cluster placed 
in the centre of the focussed beam and ii) the position of 
the cluster integrated over the approximately estimated 
spatial profile of the pulse. As we will see later, at higher 
pulse fluences averaging over the spatial pulse profile will 
have a higher impact on the model predictions, compar- 
ing to the case of a cluster placed in the centre of focussed 
beam. This is due to the strong nonlinear dependence 
of the ion charge created within a cluster on the pulse 
fluence at higher pulse fluences. At lower fluences this 
dependence is linear. 

First we will investigate ionization dynamics at differ- 
ent cluster sizes. In figs. [IJi-d we plot the charge state 
distributions at the end of expansion phase obtained for 
four different cluster sizes irradiated at a fixed pulse flu- 
ence of 0.4 J/cm 2 . The trend can be understood in the 
following way. The inverse bremsstrahlung is the domi- 
nant mechanism of energy absorption within our model. 
This process heats up the quasi-free electrons within the 
cluster. The hot electrons collide with ions and atoms, 
producing higher charges and releasing new electrons. 
These processes compete with the three-body recombi- 
nation processes that reduce the ion charge and decrease 
the electron density. After the pulse is over, and the sys- 
tem reaches the local thermodynamic equilibrium (LTE) 
state, the number of ionization and recombination events 
saturates, changing slowly with the decreasing tempera- 
ture within the cluster. 

Within small clusters a large fraction of electrons re- 
leased during the ionization processes can leave the clus- 
ters early in the exposure. The remaining electrons are 
not heated efficiently via inverse bremsstrahlung pro- 
cesses due to their low density within the cluster. Con- 
sequently, only low charge states are observed (fig.QJi-b). 
In contrast, within large clusters of 2500 and 90000 xenon 
atoms only a small fraction of the released electrons are 
able to escape from the cluster. The width of the pos- 
itively charged outer shell is small with respect to the 
radius of the neutral core. Electrons confined within the 
core are then heated efficiently. That leads to further 
collisional ionizations. In fig. QJj-d we plot the respec- 
tive ion fractions. Our predictions are in agreement with 
experimental predictions. The trend of the size depen- 
dence is correct: at fixed pulse energy the maximal ion 
charge created increases with the cluster size until it sat- 
urates at larger cluster sizes. The average charge is plot- 
ted as a function of cluster radius in fig. [T]3. Similarly, 
at higher cluster radii the average charge saturates. This 
is in agreement with the experimental data. Saturation 
of the ion charge created within large clusters irradiated 
with a pulse of a fixed fluence is due to the fact that for 
clusters large enough the energy absorbed from the pulse 
(per atom) will not be sufficient for the creation of higher 
charge states. 

Below we also show our estimates for the average ki- 
netic energy per ion, E, as a function of the cluster radius 
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FIG. 1: Ion fractions, Ni/N, within the irradiated clusters estimated for various cluster sizes: (a) 20 atoms, (b) 70 atoms, 
(c) 2500 atoms and (d) 90000 atoms. The number, N, denotes here the number of Xe +% ions, and the number, N, is the total 
number of ions. Model predictions for clusters placed in the centre of the focussed beam and predictions integrated over the 
approximately estimated gaussian spatial pulse profile are shown. Also: (e) average charge, Z, created within the irradiated 
cluster, and (f) average energy absorbed per ion, E, within the irradiated cluster as a function of the cluster radius, R. The 
model estimates (squares) were obtained with pulses of different intensities and lengths but of fixed integrated radiation flux, 
F — 0.4 J/cm 2 , and then averaged over the number of pulses. Experimental data are plotted with stars. Model predictions for 
clusters placed in the centre of the focussed beam (open squares) and predictions integrated over the approximately estimated 
spatial pulse profile (filled squares) are shown. 
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FIG. 2: Ion fractions, Ni/N, within the irradiated xenon clusters (N a toms = 2500). These clusters were irradiated with 
rectangular pulses of fluences: (a) F = 0.05 J/cm 2 , (b) F = 0.3 J/cm 2 , (c) F = 0.84 J/cm 2 , and (d) F = 1.5 J/cm 2 . The 
number, JV», denotes here the number of Xe +l ions, and the number N is the total number of ions. Model predictions for 
clusters placed in the centre of the focussed beam and predictions integrated over the approximately estimated gaussian spatial 
pulse profile are shown. Also: (e) average charge, Z, created within the irradiated xenon cluster (Natoms = 2500), and (f) 
average energy absorbed per ion, E, within the irradiated cluster as a function of the pulse fiuence, F. The model estimates 
(squares) were obtained with pulses of different intensities and lengths but of fixed fiuence, F = 0.05 — 1.5 J/cm 2 , and then 
averaged over the number of pulses. Experimental data are plotted with stars. Model predictions for cluster placed in the 
centre of the focussed beam (open squares) and predictions integrated over the approximately estimated spatial pulse profile 
(filled squares) are shown. 
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(fig- [If)- The energy absorption per atom increases with 
the cluster size and also saturates for larger clusters. 

The dependence of the ionization dynamics on the radi- 
ation flux has been experimentally investigated for xenon 
clusters consisting of 2500 atoms. Below we show the 
plots of the ion fractions obtained with the experimental 
data and the ion fractions obtained with our model (fig. 
d^-d) ■ Higher pulse fluences lead to the creation of higher 
charge states within the clusters in accordance with the 
mechanism described above. 

At the lowest radiation flux, F — 0.05 J/cm 2 (e. g. 
1 = 5- 10 12 W/cm 2 and At = 10 fs for a rectangular 
pulse), charges up to +3 were found. Our model shows 
only up to doubly charged ions. As suggested in ref. 
[221 ] . this small discrepancy between experiment and our 
predictions may be due to the fact that at low fluences the 
many-body recombination effects within the cold electron 
plasma can be important [20j ] . They could then lead to 
higher ionization states that we do not observe within 
this model. 

At the flux, F = 0.3 J/cm 2 , experimental ion fractions 
have a maximum at Xe +2 , whereas our predictions peak 
at Xe +1 . However, the maximal ion charge is found to 
be +5 with both experimental data and simulation re- 
sults. At higher flux, F = 0.84 J/cm 2 and the F = 1.5 
J/cm 2 , the ion fractions are in a good agreement with 
experimental data. Maximal charges up to +7 and +8 
are observed respectively. 

The average charge is plotted as a function of radia- 
tion flux in fig. [2^. With the increasing pulse fluence the 
average charge created increases. The charges calculated 
are very close to the corresponding experimental values. 
Also, as expected, the content of neutral atoms found 
within the cluster (not shown) decreases at the increas- 
ing pulse fluence. 

Below we show also the average kinetic energy per ion 
(estimated with our model) as a function of the radia- 
tion flux (fig. [2f)- The sparse experimental data do not 
allow the identification of the trend of the fluence depen- 
dence, i.e. whether it is linear or non-linear. Our predic- 
tions slightly underestimate the experimental predictions 
at high fluences but stay within the error limit given by 
the experimental uncertainty of fluence estimation. 

In summary, we used a microscopic model based on 
a first principle Boltzmann approach to investigate non- 
equilibrium dynamics within atomic clusters of various 
sizes irradiated by single VUV pulses. This model is 
computationally efficient for small and large clusters, and 



includes various processes that are relevant at VUV pho- 
ton energies. Predictions obtained at various cluster sizes 
and various pulse fluences were compared to the experi- 
mental data and found to be in good overall agreement, 
especially if one considers the experimental uncertainty of 
pulse fluence estimation given by a factor of 2.5. The re- 
sults obtained have been cross-checked with independent 
molecular dynamics simulations, and have also showed 
good agreement. 

The results obtained demonstrate the different ioniza- 
tion dynamics of small and large clusters. For small 
clusters the efficient escape of electrons from the clus- 
ter reduces the number of higher charges created within 
the sample during collisional ionizations. Recombination 
is also suppressed due to the absence of free electrons. 
For large clusters, many electrons stay within the sam- 
ple as they are kept there by the attractive potential of 
the outer shell. They are heated by the IB process, and 
ionize and recombine efficiently. Higher charges are cre- 
ated, and as a result of efficient recombination neutral 
atoms appear within the cluster core. After the pulse is 
over, the system reaches the LTE, and the total num- 
ber of ions and atoms within the sample changes only 
weakly. A fraction of neutral atoms remains within the 
core, its magnitude depending on the cluster size. This 
indicates that the ions detected during the experiment 
come mainly from the surface and the outer part of the 
cluster. 

In conclusion, we find that formation of high charge 
states and the strong absorption of VUV radiation is 
quantitatively understood within the framework of our 
model. 
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